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ANALYSIS OF THE THEORY OF HIGH-ENERGY ION TRANSPORT 


John W. Wilson 
Langley Research Center 

SUMMARY 

In the operation of commercial high-altitude aircraft and manned spacecraft, expo- 
sure to the heavy -ion radiations of space must be considered. There is a need, therefore, 
to understand the Interaction of these radiations with the structures and shielding mate- 
rials of these transport vehicles. The development of analytical methods to study the 
transport properties of high-energy ions in materials is the topic of the present report. 

Procedures for the approximation of the transport of high-energy ions are discussed 
on the basis of available data on ion nuclear reactions. A straightahead approximation 
appears appropriate for space applications. The assumption that the secondary -ion- 
fragment velocity is equal to that of the fragmenting nucleus is found to be inferior to 
straightahead theory but may be of sufficient accuracy if the primary ions display a broad 
energy spectrum. An iterative scheme for the solution of the inhomogeneous integral 
transport equations holds promise for practical calculation. A model calculation shows 
that multiply charged ion fragments (atomic number greater than 3) are able to penetrate 
to greater depths in comparison with the free path of a primary heavy ion (atomic number 
approximately equal to 25). 


INTRODUCTION 

The prospect of extensive space operations in the era of the Space Transportation 
System, of ever increasing altitude range of commercial aircraft operations, and of 
recent development of high-energy heavy-ion accelerators (especially in the context of 
radiotherapy) accentuates the need for a theoretical understanding of the interaction of 
energetic heavy ions with extended matter. The main limits on such developments in the 
past have been a result of the scarcity of information on the nuclear reactions Induced by 
heavy-ion collisions (ref. 1). With the advent of the acceleration of heavy ions to relati- 
vistic energies by particle accelerators (refs. 2 and 3), needed experimental data will be 
forthcoming and the further development of heavy-ion transport theory now seems appro- 
priate. In addition, progress in heavy-ion reaction theory is likewise being made (refs. 4 
to 10); many of the transport parameters may be accurately estimated and Implications of 
heavy-ion dynamics on the transport theory can now be investigated. 



It is considered that the simplest approach to studying high-energy ion transport 
would be to implement the Monte Carlo method. However, as a result of the large num- 
ber of possible outcomes of heavy-ion nuclear reactions, it is felt that inordinately large 
samples will be required to adequately represent even the more important quantities of 
the transport process. For this reason, the development of analytical methods appears 
to be a potentially fruitful area of research. 

The purpose of the present paper is to make a preliminary analysis of heavy-ion 
transport theory with particular attention to the development of analytical methods well 
suited to numerical approximation without recourse to Monte Carlo methods (ref. 11). 

Since radiation protection is of concern, consideration will be given to evaluation of factors 
closely related to biological response (ref. 12). 

SYMBOLS 

Aj atomic mass of type j ion, amu 

E ion kinetic energy, MeV/amu 

fjjj(E,E',f2,fi') probability density of type k projectile with energy E' and 

motion along to produce type j ion with energy in dE 
about E and in dn about fi, (sr-MeV/amu)"^ 

G-(x,fi,E) source of type j ions at x with energy E and direction fi, 

J 4 

(sr-cm^-MeV/amu)' 
m nucleon mass, amu 

if(r ) outward directed imit normal vector at point F on bounding surface 

Pj(E) total survival probability of type j ion of energy E 

Rj(E) continuous slowing-down range of type j ion of energy E, cm 

rjj^(E,E',fi,fi') double differential cross section for production of type j ions of 

energy E into direction by type k ions of energy E’ from 
direction (cm-sr- MeV/amu)"^ 

Sj(E) stopping power or linear energy transfer (LET) due to interaction of type j 

ion with orbital electrons of transport medium, MeV/cm 
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s,t two-parameter net defined over bounding surface 

X one -dimensional position vector, g/cm^ 

X position vector (see fig. 1), cm 

Xj! component of x perpendicular to fi, cm 

position vector defining point on boundary which is origin of vector connecting 
boundary to x along direction n (see fig. 1), cm 

y(s,t) two-parameter vector function defining the bounding surface, cm 

i^jj^(E) multiplicity of type j ions produced by collision of type k ion of energy E 

|j,T 7 j characteristic coordinates of type j ion, cm 

p projection of x onto n, cm 

Sj(E) macroscopic absorption cross section for type j ion of energy E, cm"^ 

CTjjj standard deviation of momentum of type j ions produced in collision of 

type k ions, MeV/c 

0j(x,S2,e) differential flux of type j ions at x directed toward with 
energy E, (cm^-sec-sr-MeV/amu)"^ 

n ion direction of motion 

Superscripts; 

P projectile 

T target 

Primes indicate a variable of summation or integration. 
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TRANSPORT THEORY 


A heavy ion, after entering a region filled with ordinary matter, interacts with 
orbital electrons, causing ionization and excitation of the medium. Because of the large 
mass difference between the ion and these orbital electrons, only a small amount of the 
ion energy can be transferred in a collision with a single electron. As a result of the 
long range of the Coulomb force and the large percentage of the material volume occupied 
by electrons, the electron interactions can, to a good approximation, be treated as a con- 
tinuous slowing -down process over any finite path lei^th. Althoi^h the energy lost by an 
ion over some fixed path length fluctuates about a mean value, this fluctuation amounts to 
no more than a few percent (refs. 13 and 14) and is of no importance in the study of space 
radiation (ref. 15). In the following, continuous slowing-down theory will be assumed 
throughout and the relevant quantity is the average energy loss per unit path length, 
denoted by Sj(E), where E is the ion energy and j denotes the ion type. 

The mean free path for nuclear collisions is large (more than a centimeter); by com- 
parison, the mean free path for collision of the ion with electrons is small. Whereas 
collisions with electrons result only in a small transfer of energy compared with the total 
ion kinetic energy, the nuclear collision generally alters (loss of mass and charge) the ion 
and the struck nucleus, with many secondary particles being produced. The secondary 
particles produced as fragments of the primary heavy ion will have longer ranges and 
free paths causing much greater penetration. As the secondaries undergo additional 
nuclear reactions, more secondaries, which penetrate deeper into the material, are pro- 
duced. The purpose here is to develop the theoretical understanding of the transport of 
such radiations in extended materials. 

The heavy-ion transport equations are derived by balancing the change in the ion 
flux as it crosses a small volume of material with the gains and losses due to nuclear 
collision (ref. 11). The resulting equations for a homogeneous material are given by 

= ^ y dE' dj!’ r.j^(E,E’,fi,S5’) <t>^(x,fi\E') (1) 

k 

where 0j{xjfi,E) is the flux of ions of type j with atomic mass Aj at x with 
motion along n and energy E in units of MeV/amu, Sj(E) is the corresponding 
macroscopic cross section, Sj(E) is the linear energy transfer (LET), and 
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rjj^(E,E’,fi,0') is the production cross section for type j ions with energy E and 
direction by the collision of a type k ion of energy E' and direction The 

term on the left side of equation (1) containing Sj(E) is a result of the continuous 
slowing -down approximation, while the remaining terms of equation (1) are seen to be the 
usual Boltzmann terms. The solutions to equation (1) exist and are unique in any convex 
region for which the inbound flux of each ion type is specified everywhere on the bovmding 
surface. If the boundary is given as the loci of the two-parameter vector function y(s,t) 
for which a generic point on the boundary is given by F, then the boundary condition is 
specified by requiring the solution of equation (1) to satisfy 

0 .(f,n,E) =)//.( f,n,E) (2) 

for each value of S2 such that 



(3) 


where n(r) is the outward directed unit normal vector to the boundary surface at the 
point r and i//j is a specified boundary function. 

Equation (1), with the boundary condition (eq. (2)), may be written as an inhomogen- 
eous integral equation. To do so requires inversion of the differential operator contained 
in brackets on the left-hand side of equation (1). As will be shown, there is a character- 
istic space in which the inversion is simplified. The transformation of equation (1) to 
this characteristic space will now be made. (See ref. 11 for a more detailed discussion.) 
Let the secondary particle source term on the right-hand side of equation (1) be denoted 
by Gj(x,n,E), so that 


■ r ^ S.(E) + 2j(E)j 0. (x,n,E) = G.(x,0,e) 


Multiplying equation (4) from the left by Sj(E) results in 


S.(E) a 

^ J ^ 


0j(x,n,E) =G.(%n,E) 


where 


0j(x,n,E) =Sj(E) 


(4) 


(5) 


( 6 ) 
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and 


g^(x,s!,e) =Sj(e) g^(x,j5,e) 


iV 


A set of characteristic variables are found to be 


p = n • X (8) 

Xn = X - pn (9) 

= P - Rj(E) (10) 

^j=P + Rj(E) (11) 

and the new functions in this characteristic space are defined as 

= + (12) 

~ Gj^pO + “SCj^jOjE^ ( 13 ) 

= 2j(E) (14) 


When equation (5) is mapped into the space defined by the variables given by equa- 
tions (8) to (11) and the new functions given by equations (12) to (14), one obtains the sim- 
plified differential given by 




(15) 


which can be formally solved. Equation (15) may be integrated using the appropriate 
integrating factor, resulting in 


<j(,j.{j,x„,s5) . e==p[-| Jj' vfi) 
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where a is any real number. If the flux ^that is, known at some 

point rj. given by the set of values r]. = a, and fi, then the flux at any other rj. 

J J J J 

may be found by using equation (16). The value a will be chosen later to always corre- 
spond to the boundary y(s,t). Equation (16) may now be rewritten with the aid of equa- 
tions (8) to (14) as 


(p.(x,n,E) = exp 


p-Rj(E) 


R. 

J \ 2 



J — n + x„,n,R. 

2 ’ J \ 2 


n 


p-Rj(E) 


1 

exp^ - \ 

2 *^p-Rj(E) J 


“i -V- 


On' 


xG. 


J + Xn,fi,Rj 


d?]" 


(17) 


The value a will now be chosen such that the second factor of the first term on 
the right-hand side of equation (17) corresponds to the value at the boundary. Hence, 
a is chosen such that 


^3 ^ ^ + X - r 

2 + ^n - ^n,x 


(18) 


The boundary point is shown in figure 1. The corresponding value of a is given by 


a - 2n • - p - Rj(E) 


(19) 


Using in equation (17) the variable transformations given by 


E’ 




( 20 ) 


and 


_l|p + Ri(E) - J?"! 

= ^3 1—^1 J 


( 21 ) 


7 


so that 


S.(E') 

dE' = 

2Aj 


( 22 ) 


and 


dE" = - 


S.(E") 

J drj" 

2Aj 


(23) 


results in 


<^.(x,f2,E) = exp 


Aj2j(E») dE' 


I, 


-f 

R.'^[Rj(E)-d+#0 


Sj(E’) 




exp 


pE" A.S.(E') dE' 
Je Sj(E') 


X Gj|(p + Rj(E) - Rj(E"i)n + \,fi,E"J^ 


A, dE" 


(E") 


(24) 


where 


d = • r. 


n,x 


Defining the total nuclear survival probability as 


P^(E) = exp 


pE S.(E') dE' 

3 Jo Sj(E') 


(25) 


(26) 


results in 
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( 27 ) 


<.j(x,a,E) = -i i <., (f„_j,n,Rj'’p - d * Rj(Egj 

] J 


Equation (27) can be further reduced using equations (6) and (7) to 







Sj(E)Pj(E) 


Gj|jc + [R.(E) - R.(E’gj!,fi,E’^ dE’ (28) 


where 


E^=R.-lf.d.Rj(Eg 

Using the definition of Gj (x,f?,E) (eq. (4)) in equation (28) results in 




Y' '^Ev A.P.(E') ^0° p ^ / — — \ 

+ ) ^ dE' J J . ■ \ dE” \ dJ2'r., (E’,E”,n,n’) 

L Jv S.tElP. El Jf* J jkV ' ’ ’ ' 


^ Sj(E)Pj(E) Je 


(29) 


X 4>. 


IjE +0t.(E) - R.(E’gj7,j2’,E"j 


(30) 


which is the integral equation of heavy -ion transport theory. Techniques for the approxi- 
mate solution of equation (30), or equivalently equation (1), with the boundary condition 
(eq. (2)) are considered next. 
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INTERACTION PARAMETERS 


The parameter appearing in equation (1) relating to interaction of ions with orbital 
electrons is the energy loss per unit path length, denoted by Sj(E), For the present, 
fluctuations about the mean will be ignored and the simplified approximation (see refs. 16 
to 18 for discussion) 


47re^z.^ 

Sj(E) = 1- N„B 

will be used, where 


B = z 


m 


/2m v^\ / iz.e^^ 


(31) 


(32) 


i//(x) = 


1 d 
r(x) dx 


r(x) 


(33) 


Njjj is the atom density of the material, Zj^^ is the material atomic number, I is the 
material adjusted ionization potential, m^ is the electron mass, e is the electron 
charge, h is Planck's constant, v is the ion velocity, and zj is the ion atomic num- 
ber. Equations (31) to (33) apply at large ion velocities. As the ion slows, it tends to 
capture orbital electrons, which effectively reduces the ion charge. The effective ion 
charge Zj will be taken as (ref. 16) 





1 - exp 1-125 v/cz 


2/3 


J 


(34) 


where v is the ion speed and c is the speed of light. In addition to orbital electron 
capture, the nuclear elastic scattering also becomes important at very low energies. If 
the ion speed is svifficiently large, it is clear from equations (31) through (33) that 


Sj(E) 



Sp(E) 


(35) 
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where Sp(E) is the proton stopping power and Zp is the effective proton charge. 
Empirically, it is found that equation (35) is reasonably accurate even at low ion speeds. 
Expressions (31) to (33) fail at low energies but the empirical expression (for energy loss 
in water) 


Sp(E) » e0*3°3(2517 - 6283E) (36) 

and equation (35) may be used for energy E less than 243.8 keV/amu. 

The parameters in equations (1) and (30) which are least known are those associated 
with the nuclear attenuation and fragmentation. If the elastic scattering of heavy ions is 
ignored to a first approximation, then Sj(E) is the ion macroscopic absorption cross 
section for the transport medium. The microscopic cross sections from which Sj(E) is 
determined are presented in reference 10. 

The fragmentation of the ion and target nuclei is represented by the quantities 
rjj^(E,E',J2,n') which are composed of three functions as follows: 

rjk(E,E’,n,n’) = Sj^(E')i^jij(E’)f.j^(E,E’,n,n') (37) 


where ^'jij(E') is the multiplicity of type j ions being produced by a collision of a 
type k ion of energy E’ and f^k(E,E’,n,n) is the probability density distribution for 
producing type j ions of energy E into direction SI from the collision of a type k 
ion with energy E' moving in direction S2'. For an unpolarized source of ions and 
unpolarized targets, the energy-angle distribution of reaction products may be taken as a 
function of the energies and cosine of the production angle with respect to the incident 
beam direction. The secondary multiplicities i^jij(E) and secondary energy-angle dis- 
tributions are the major unknowns in ion transport theory. 


Until recently, information on the multiplicity i'jjj(E) was obtained through exper- 
iments with galactic cosmic rays as an ion source and the fragmentation of the ions on 
target nuclei was observed in nuclear emulsion (ref. 19). Such data are mainly limited 
by not knowing the identity of the initial or secondary ions precisely and by relatively low 
counting rates of each ion type. Recent heavy-ion acceleration by machine makes it pos- 
sible to reduce the uncertainty since large count rates are possible with known ion types. 
In addition, the target nuclei in accelerator experiments can conveniently be other than 
nuclear emulsion, and accurate detector techniques with modern electronic processing are 
greatly improving the experimental data base. In addition, the accelerator experiments 
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are providing information on the spectral distribution which has not before 

been available (ref. 20). 

The spectral distribution function is found to consist of two terms which describe 
the fragmentation of the projectile and the fragmentation of the struck nucleus as follows 
(refs. 21 and 22): 


rjj^(E,E',j2,fi’) = Sj^(E') (E')fj^(E,E’,n,fi’) + i.J(E’)f J(E,E’,n,fi'^ (38) 


P P T T 

where and fj^ depend only weakly on the target and and fj^^ depend only 

p 

weakly on the projectile. Whereas the average secondary velocities associated with f 

T 

are nearly equal to the projectile velocity, the average velocities associated with f are 
near zero. Experimentally (refs. 20 to 22), it is observed that 


fj^(E,E',i!,f2-) * 


m 


"("ik)' 


3/2 


/2E 


exp 


. (p - p') ^ 


2 a.' 


jk 


m 


2J7^a.,.\2 


^Jk)' 


3/2 


\f^ exp 


(/2mEfi - \/2mE' 
^2 


2'-3^k 


(39) 


(where p and p' are the momenta per unit mass of j and k ions, respectively) 
and 


f-^(E,E’,fi,fi’) - 


m 




3/2 


^/2¥ 


exp 


P^ 


2 a. 


jk 


(40) 


P T 

where ajj^ and aj^ are related to the rms momentum spread of secondary products. 

These parameters depend only on the fragmenting nucleus. It was first suggested by 
Feshbach and Huang (ref. 7) that the parameters a^ and aj^ depend on the average 

square momentum of the nuclear fragments as allowed by Fermi motion. A precise form- 
ulation of these ideas in terms of a statistical model was obtained by Goldhaber (ref. 8). 
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APPROXIMATION PROCEDURES 


Neglect of Target Fragmentation 


Employing equations (37) to (40) in the evaluation of 
results in 



of equation (4) 


Gj(x,i2,E) = ^ JdE' dJ2' Sj^(E’) 0j^(x,n’,E’) 
k 

= g?(x,o,e) +G-^(x,n,E) 





(41) 


where, as before, the superscripts P and T refer to fragmentation of the projectile 
and target, respectively. The target term is seen to be 


gT(x,J2,e) 


m 




3/2 


\[m 


exp 


J dO’ dE’ i/^(E’) 2k(E’)0j^(x,s!’,E') 

(42) 


which is negligibly small for 

^2/ 


E » (cr.j^ 1 /m 


(43) 


Thus, for calculating the flux at high energy, 



(44) 


Space Radiations 


A convenient property of space radiations is that they are nearly isotropic. This 
fact, coupled with the forward peaked spectral distribution, leads to substantial reductions 
in the source term as follows: 


G?(x,n,E) JdE’ dsT’ Sjj(E') i/^(E') 


m 


l3/2 






X exp 


(\/2mE?l - \/2mE'n*) 


0j^(x,fi’,E’) 


(45) 
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Assuming that is a slowly varying function of fi', one may seek an expan- 

sion about the sharply peaked maximum of the exponential function. Such an expansion is 
made by letting 

O' = S2 + (cos 6 - 1)S2 + sin 6 "e^ (46) 

where 


cos 0 = O • 


and 


e,=nx O'y 

<p 


with which the flux may be expanded as 


0k(x,fi’,E') = 0j^(x,n,E') + 


^ 0j^(x,n,E')' 


9fi 


(cos 6 - 1)S2 + sin 6 "5 


4 > 


+ 


(47) 


(48) 


(49) 


Substituting equation (49) into equation (45) and simplifying results in 



The leading term of equation (50) is clearly a good approximation to the source term 
whenever 


2mE 



» 


dQ, ^ 


(51) 


Note that the leading term is equivalent to assuming that secondary ions are produced 
only in the direction of motion of the primary ions. In the case of space radiations which 
are nearly isotropic, relation (51) is easily satisfied and neglect of higher-order terms 
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in equation (50) results in the usual straightahead approximation. If the radiation is 
highly anisotropic, then relation (51) is not likely to apply. Such behavior was discovered 
empirically by Alsmiller et al. (refs. 23 and 24) for the case of proton transport. 

Velocity Conserving Interaction 

It has been customary in cosmic ion transport studies (ref. 1) to assume that the 
fragment velocities are eqiaal to the fragmenting ion velocity prior to collision. Derived 
below is the order of approximation resulting from such an assumption. If it is assumed 
that the projectile energy E* is equal to the secondary energy plus a positive quantity e 

E' = E + e (52) 


and that e will contribute to equation (50) only over a small range above zero energy, 
substituting equation (52) into equation (50) and expanding the integrand results in 



Since \ « 1 at those energies for which most nuclear reactions occur, it is 

V mE 

clear that the assumption of velocity conservation is inferior to a straightahead approxi- 
mation but may be adequate for space radiations where the variation of 0j^(x,n,E) with 
energy is sufficiently smooth. That is, 

^ ^ ~ <|)k(x,n,E) 


Decoupling of Target and Projectile Flux 
Equation (1) with the use of equation (41) may be rewritten as 

B.0.(x,fi,E) F^ 
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where the differential operator is given by 


Bj = 


. V - i-A S.(E) + s. 
Aj 9E ' 


j(E) 


and the integral operator (^jk - ^^k + is given by 

F.^<t>y.{x,fi,E) = J dE’ dfi'r.j^(E,E’,n,s5’) 0j^(x,fi’,E’) 


(55) 


(56) 


Defining the flux as a sum of two terms 


0.(x,fi,E) = 4>J{x,Q,e) + <^?(x,0,e) 


(57) 


allows the followii^ separation: 


B.(pf(x,a,E) =J^Ffy.<l>l{xn,E) Ff^(f>l{%n,E) 


(58) 


B.cpJ {%n,E) Fj0j(x,fi,E) 


(59) 


As noted in connection with equations (42) through (44), the source term on the right-hand 
side of equation (59) is small at high energies and one may assume 


«o 


(60) 


for E » ^m. As a result of equation (60) and the fact that the ion rai^e is small 

compared with its mean free path at low energy, one obtains 


Bj(^?(x,fi,E) *^Fj^0P(x,J2,e) 


( 61 ) 


B.4,J(x,a,E) * J fJ0^(1c,S2,e) 


( 62 ) 
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The advantage of this separation is that once equation (61) is solved by whatever means 
necessary, then equation (62) can be solved in closed form. The solution of equation (62) 
is accomplished by noting that the inwardly directed flux 4>J must vanish on the bound- 
ary, so that 



AjPj(E') 

Pj(E)Sj(E) 


y dE" dfi'rj^(E’,E",jl,n’) 


X 


+ (R.(E) - Rj(E 



(63) 


where E^ is given by equation (29). 

Using equations (38) and (40) in equation (63) yields 



AjPj(E') 


Pj(E)Sj(E) 


m 


2jt 




3/2 





where 



+ [Rj(E) 


Rj(E 



?j(x) = ^ J clE’ dS2'Sj^(E') i^J^(E') 0^(x,fi',E') 
k 


(64) 

(65) 


T 

and has been assumed to be a slowly varying function of projectile type k and pro- 

jectile energy E. K the range of secondary type j ions is small compared with their 
mean free path lengths and the mean free paths of the fragmenting parent ions fj^, that is, 

« fk (66) 



then the integral of equation (64) may be simplified as 



(67) 
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which may be reduced into terms of known fvinctions. Thus, 



in terms of the incomplete gamma function. It may also be shown that equation (68) is 
equivalent to 



At points sufficiently removed from the bovindary such that 



equation (69) may be reduced to 



The solution of equation (61) will now be further examined. 

Back-Substitution and Perturbation Theory 

One approach to the solution of equation (61) results from the fact that the multiply 
charged ions tend to be destroyed in nuclear reactions. Thus, 

=0 (j s k) (72) 

This means that there is a maximum j such that 
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I 


Bj0j(x,n,E) = 0 


(73) 


where J is the largest j. Furthermore, 


(xjOjE) 


(74) 


and, in general. 


N-1 

k=l 

for N < J - 1. Note that equations (73) to (75) constitute solvable problems. The sii^ly 
charged ions satisfy 

J 

Bi0f(x,n,E) = F^^J’(lc,j2,E) + ^ F^j^(/,f(x,?2,E) (76) 

k=2 


which, unlike equations (73) to (75), is an integral-differential equation, which is difficult 
to solve directly. Equation (76) is solvable by perturbation theory and the resultant 
series is known to converge rapidly for intermediate and low energies (ref. 11). It is to 
be noted that equations (73) to (75) are also obtained from perturbation theory as applied 
to equation (61) at the outset. Thus, the perturbation series is expected to converge 
after the first J plus a few terms. 

When J is large, the number of terms enterii^ the summations in equations (75) 
and (76) may lead to a lengthy computational problem. An estimate of the number of 
terms required to adequately approximate the solution of equation (75) would be useful 
before a commitment to the development of computer programs to evaluate the series is 
made. For this purpose, a study of a simplified but realistic model is undertaken to 
determine the rate of convergence of the perturbation series. 

MODEL CALCULATION 

At sufficiently high energies, the energy loss can be neglected over many mean free 
paths of the incident ion. Assuming the straightahead approximation is valid and that all 
secondary fragments have the same velocity as the fragmenting parent nucleus leads to 
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a particularly simple one -dimensional form of the theory. Namely, equation (1) becomes 




which may likewise be written as 


(77) 


-S.X pX p 

’#')(*) = * •t’P * Jo ® Z, ^k‘'jk'»'k(=‘'> *=' 


X S.(x*-x)' 


(78) 


which is equivalent to equation (30). For the present example, a simple form for 


and which is reasonably realistic is assumed. Namely, we take 

2k = 2^/8 


(79) 


and 


where j and k are to be interpreted as ion charge numbers and u is the unit step 
function; it is seen that equation (80) conserves charge in the interaction. Implicit in 
equations (79) and (80) is the neglect of target fragmentation products. 

Equation (78) was solved numerically using successive interations and numerical 
quadrature. The zeroth-order approximation was taken as 


^j°\x) = e ^ 0.(0) 


(81) 


and the higher-order corrections were added by substitution into the right-hand side of 
equation (78), giving 


0S'-"^)(x) = e 


-S,x 


03-(O) 


■2 


p px S.(x’-x) -X 


(82) 


for i = 0, 1, 2, . . . . The 2^ parameter in equation (79) was chosen to correspond to 
air and the boundary condition was taken as 
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0^(0) = 25 


(83) 


corresponding to an iron group nucleus incident on the Earth's atmosphere. The iteration 
procedure converged over the first 60 g/cm2 in six Iterations and numerical errors were 
held to within 10 percent. Thus, it appears that perturbation theory may well result in a 
practical calculation procedure for heavy-ion transport for most applications. 

The integrated ion flux as calculated with the iteration procedure indicated in equa- 
tion (82) is shown in figure 2, along with the integrated ion flux as calculated with only 
the first-order correction (that is, i = 0 in eq. (82) only). It is seen that the first-order 
result is far from accurate over most of the regions of interest. Although first-order 
theory is inadequate in this respect, the predicted average charge of the beam as it passes 
through the air is no worse than a factor of two off, as shown in figure 3. Although the 
average charge is approximately correct, 96 percent of the initial charge of 25 is con- 
tained in higher-order terms at 60 g/cm2. The importance of higher-order terms is 
shown again in connection with the calculation of secondary particle dose in figure 4. 
Clearly, doses calculated on the basis of first-order theory contain large errors for 
depths beyond several g/cm2 of material. 

The dose contribution from various charge components is shown in figime 5 at sev- 
eral depths. Also shown for comparison are the results of first-order theory for the two 
lowest values of depth. Most conspicuous is that the maximum contribution comes from 
a broad range of charges above z = 3. This is even true at rather great depths such as 
the 59.2 g/cm2 curve in the figure, which represents nearly four mean free path lengths 
of the primary ion beam. 

The results in figure 5 are of importance in the evaluation of the biological effects 
since it is well known that heavy ions can bring about biological effects not observed with 
any other known types of radiation. The fact that secondary ions formed from the frag- 
mentation of a single larger ion have greater penetrating power and may be even more 
damaging biologically as a result of saturation effects points out the need for a careful 
evaluation of transport effects on biological systems shielded from high-energy heavy -ion 
radiation. Such an evaluation would also be important with respect to evaluation of doses 
in connection with radiotherapy (especially in surrounding healthy tissues). 

Aside from the damaging effects of secondary projectile fragments, the nuclear 
stars formed by fragmenting target nuclei will sometimes be of major importance. In 
this connection, the collision density is of importance as shown in figure 6. Again it 
is seen that a first-order calculation does not produce the required accuracy for most 
applications. 
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CONCLUDING REMARKS 


The fvindamental equations of heavy -ion transport are presented both as a set of 
integral/partial-differential equations and as a coupled set of inhomogeneous integral 
equations. The basic interaction parameters are discussed and simple approximations 
based on experimental observations are given. On the basis of these interaction param- 
eters, simplifications of the fundamental transport equations are discussed. An iterative 
method is proposed to solve the inhomogeneous integral equations for which rapid conver- 
gence is obtained in a transport model which neglects ionization energy loss, assumes 
forward secondary particle production, and assumes the secondary particle velocity is 
equal to the fragmenting ions velocity. Although a first-order perturbation theory is 
generally insufficient to represent the transport process, six correction terms provide 
svifficient accuracy for the iron group ions for most depth ranges of interest (^60 g/cm^). 
It is found in these model calculations that the secondary multiply charged ions can pene- 
trate well beyond the depths obtainable by the primary ion. This fact will be most impor- 
tant in connection with space and high-altitude aircraft shielding and radiotherapy. 

Langley Research Center 

National Aeronautics and Space Administration 
Hampton, VA 23665 
November 24, 1976 
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Integrated ion flux, ions/cm 



Figure 1.- Relation of boundary and vector quantities required 
for evaluation of ion fields. 



X. g/cm^ 

Figure 2.- Integrated ion flux as result of normally incident iron group 
ion flux on top of Earth's atmosphere. 





Figure 3.- Average charge of ions produced in Earth’s atmosphere by 
normally incident high-energy iron group ion integrated flux. 



X, <3 /cm 

Figure 4.- Dose in Earth's atmosphere due to normally incident iron group 
integrated flux (l/cm^) of high energy. 




Figure 5.- Dose contribution from various ion types in Earth’s atmosphere 
due to normally incident high-energy iron group ion integrated flux (l/cm2). 



Figure 6.- Density of nuclear reactions with air nuclei produced by normally 
incident iron group ion integrated flux (l/cm2) and all secondary fragments. 
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